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Abstract: Recently it was shown that there is a class of models in which colored vector and 
scalar resonances can be copiously produced at the Tevatron with decays to multijet final 
states, consistent with all experimental constraints and having strong discovery potential. 
We investigate the collider phenomenology of TeV scale colored resonances at the LHC and 
demonstrate a strong discovery potential for the scalars with early data as well as the vectors 
with additional statistics. We argue that the signal can be self-calibrating and using this 
fact we propose a search strategy which we show to be robust to systematic errors typically 
expected from Monte Carlo background estimates. We model the resonances with a phe- 
nomenological Lagrangian that describes them as bound states of colored vectorlike fermions 
due to new confining gauge interactions. However, the phenomenological Lagrangian treat- 
ment is quite general and can represent other scenarios of microscopic physics as well. 
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1. Introduction 



Hadron colliders provide the highest center of mass energy in particle physics and are our 
portal to discovering heavy colored particles. We are currently at the dawn of the LHC era 
where we hope to discover not only the Higgs boson but other states beyond the Standard 
Model (SM) as well. While the "hierarchy problem" has been a useful guide in approaching 
this critical time, one should be open to the possibility that the mechanism responsible for the 
hierarchy may be beyond the reach of the LHC, and the data may present us with surprises 
that may or may not be related to this mechanism. From a more general viewpoint, it is 
important to look for other principles to direct us towards the discovery of new physics in 
order to make sure that we can take full advantage of the discovery potential of the LHC. 

A general strategy in approaching the unknown is simply to utilize the main strength of 
one's tools, and in the case of hadron colliders this means the production of colored states 
with large cross sections. With this in mind, it is surprising that there have been no dedicated 
searches for new physics in the multijet channel in the past. This is even true for models of 
new physics that do have multijet final states, but since the overwhelming majority of such 
models are aimed at solving the hierarchy problem they almost always contain additional 
final states that can be used to reduce background, such as leptons or the presence of missing 
transverse energy. The obvious difficulty for a search in the multijet channel is the intim- 
idating magnitude of the QCD background, however it has recently been shown in ref. [1] 
that there is a broad class of models beyond the SM with minimal ingredients, where the 
new physics would manifest itself only in the multijet channel but nevertheless the signal can 
easily compete with the background in terms of cross section. These models only require the 
presence of new colored fermions (which already appear in almost all proposed extensions of 
the SM and are essential for their discovery at hadron colliders) that are also charged under 
an additional gauge group which confines at energies moderately large compared to the QCD 
scale. The bound states of the new gauge group, dubbed "hypercolor" , can then be produced 
copiously at hadron colliders, their existence is consistent with all experimental bounds even 
if they have masses as light as a few hundred GeV, and most importantly they can be discov- 
ered with a dedicated resonance search in multijets. It was shown in ref. [1] that the Tevatron 
has a strong discovery potential for such models at sub-TeV scales, and the aim of this paper 
is to provide a complimentary outlook for an LHC search in the TeV mass range. 

When colored new fermions get bound into mesons of hypercolor, the state that cor- 
responds to the rho meson of QCD is a color octet vector particle called the coloron.^ The 
coloron can be resonantly produced at hadron colliders as it can mix with the SM gluon. This 
process can microscopically be viewed as the pair creation of the new colored fermions which 
get bound by the strong hypercolor interactions. In ref. [1] it was shown that if hypercolor 
is taken to be SU(3) with a set of fermions that transform as fundamentals of SM color as 
well as hypercolor, then this model becomes a clone of QCD in terms of its ingredients, and 
one can use QCD as an analog computer to sidestep the computational limitations of strong 

^We borrow here the terminology of refs. [2] and [3]. 
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dynamics and quantitatively predict the spectrum of the hypercolor mesons. In fact, for this 
benchmark model as well as for many other choices of the hypercolor gauge group and the 
representations of the matter fields, the coloron decays to a pair of color octet pseudoscalar 
mesons of hypercolor, called hyperpions, which in turn decay to a pair of gluons, leading 
to events with a four jet resonance. These states also appear in non-minimal Technicolor 
theories, sec for example refs. [4, 5, 6]. 

In this paper we will begin by retracing the analysis of ref. [1] for the LHC, where the 
pair production of hyperpions through their minimal QCD gauge coupling becomes more sig- 
nificant than their resonant production through the coloron, due to the dominance of gluon 
PDF's at the LHC over quark-antiquark PDF's. While the strong evidence for the secondary 
resonances is unaffected, this can weaken the evidence for the presence of the primary reso- 
nance, the coloron. Therefore we will extend the analysis of collider phenomenology for such 
models by including the process of coloron pair production, which leads to eight jet final 
states. Even though current Monte Carlo based methods can be unreliable in estimating de- 
tailed kinematic features of high multiplicity multijet backgrounds, we will argue that signal, 
where the primary production process has a two body final state, has a large enough cross 
section to overcome background involving the QCD production of hard and well separated 
multi-parton final states, even without taking advantage of the detailed kinematic features. 
In addition, we will demonstrate that signal events can be reconstructed to reveal kinematic 
features which have very low probability to be present in background events, thus leading 
to a further very significant enhancement of signal over background and making the coloron 
discoverable. For the sake of completeness we will report the results of applying our recon- 
struction procedure on a background sample that is plausibly realistic, thereby providing a 
check on our conclusions. 

The clear invariant mass peak in hyperpion pair production and the large signal cross 
section for coloron pair production also allows for a simple search strategy where the signal 
can be self calibrating. We will propose a guideline for such a search where the coloron mass 
scale can manifest itself as an anomaly in the hx distribution after requiring eight hard jets. 
We will argue that this procedure can work even allowing for relatively large systematic errors 
on the background. This measurement can be used to optimize the cuts for the subsequent 
part of the analysis for coloron mass reconstruction. 

The paper is organized as follows: In section 2 we will outline the general scenarios of 
new physics we are interested in that lead to multijet resonances at the LHC, and specify a 
phenomenological Lagrangian that describes the production and decay of the hypermesons. 
The parameters of this phenomenological Lagrangian are fundamentally determined by the 
choice of the hypercolor interactions and the transformation properties of the new matter 
fields under SM color and hypercolor. We will then concentrate on a representative choice 
for these parameters and use this benchmark model to show in section 3 that the LHC has 
strong discovery potential for the scenarios at hand, with a relatively small amount of data 
for a wide range of mass scales. We will demonstrate that combined evidence from four jet 
and eight jet final states can conclusively point to the new physics scenarios described above. 



- 2 - 



Finally, we will outline a search strategy where signal can be self calibrating such that each 
stage of the analysis determines the optimized cuts to be used at the next stage. We will 
present our conclusions in section 4. 

2. A Quantitative Model 

The scenarios we wish to discuss require a minimal set of ingredients in addition to the SM, 
the Lagrangian can be written down as 

C = CsM + i'{iIp-m)ij-^-H^,H'''' . (2.1) 

where the ^ are fermions charged under SU{i)c as well as the hypercolor gauge group (rep- 
resented by Hnv). Despite the simplicity of this setup, the resulting phenomenology is rich, 
and was discussed in detail in ref. [1]. While the choices for the hypercolor group and the 
group representations that the ijj transform under are manifold, certain qualitative features 
of the spectrum below Ahc^ the hypercolor confinement scale are common. For instance, the 
existence of an SU{3)c adjoint vector p^^ (the coloron) and the parametrically lighter SU{3)c 
adjoint scalar tt (hyperpion) are robust (we use tildes throughout the paper to distinguish 
bound states of hypercolor from those of QCD). Even though we can say little else in quantita- 
tive detail about many of the possible choices, specifying the tp to be SU{3)c triplet fermions, 
these states are the only ones below Afjc relevant for collider phenomenology. Below, we 
write down a phcnomenological Lagrangian that captures the interactions of these states as 
well as their couplings to the SM, expanding upon the one presented in ref. [1] by adding two 
terms relevant for the LHC phenomenology that were unimportant for the Tevatron. ^: 

^es = + qir {d^ + i93 (G^ + ep^)) q 

1 

+iXgsTr {G,, [p^, p^]) + ^?!^£^^ IV {{D'^p, - D,~^) [G/, G/] ) (2.2) 

m- 

The first line contains the SM kinetic terms as well as the coupling of the coloron to 
the quarks, parameterized by the effective mixing parameter e. The second line contains the 
quadratic coloron terms, where stands for the S'J7(3)c covariant derivative acting on an 

adjoint. Note that both the coloron and hyperpions have minimal gauge couplings to gluons 
due to their color adjoint nature, this allows them to be pair produced from a g-g initial 
state, which we will refer to as "pair production through QCD" . The terms on the third line 

^Wc assume here that the coloron and gluon kinetic terms have already been diagonalized, for further 
details see ref. [1]. 
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were not included in ref. [1] as they were not relevant for the Tevatron phenomenology, they 
are important for the LHC however and we include them here (with the SU{3) generators 
normalized as Tr (T^T^) = These terms represent strong interaction matrix elements 

of the underlying theory which cannot be extracted from the QCD analogy, we do how- 
ever write them down making their natural sizes manifest in the spirit of naive dimensional 
analysis (NDA) [7, 8] with free coefficients of order one. The first term of the third line is 
renormalizable and is gauge invariant by itself, with a coefficient we will say more about 
this shortly. The second term of the third line is non-renormalizable, however it is the leading 
term that leads to resonant coloron production from two incoming gluons, and since the q-q 
PDF's are not as dominant at the LHC as at the Tevatron, wc choose to keep this term. We 
estimate the size of its coefficient using NDA up to an 0(1) number, which we denote by 
The last line includes the kinetic terms for the hyperpion, the effective vertex that leads to 
the dominant coloron decay mode into a pair of hyperpions parameterized by gp-fj-jj- as well as 
the axial anomaly term that leads to the hyperpion decay into two gluons. 

At this point the effective parameters appearing in eq. (2.2) are undetermined. Although 
some of them are in principle calculable on the lattice, only little is known for generic choices 
of the hypercolor group or the hypercolor representation of the tjj.'^ There is a special case 
however, where almost all parameters can in fact be determined quantitatively. This special 
choice corresponds to the SU{3) hypercolor gauge group, with the ip transforming in the 
fundamental representation. Then the zeroth order hypercolor dynamics below Ahc (treating 
QCD as a perturbation on hypercolor) are identical to the zeroth order QCD dynamics below 
Aqcd (treating QED and quark masses as a perturbation on QCD). In the remainder of the 
paper we will refer to this choice as the benchmark model. Ref. [1] extracted the effective 
parameters of the phenomenological Lagrangian in eq. (2.2) in the benchmark model by 
scaling up the QCD spectrum and the effective couphngs of the QCD chiral Lagrangian to 
obtain 

£-0.2, (2.3) 

9m ^ 6 , (2.4) 

^-0.3, (2.5) 
nip 

/t (2.6) 



A-HC Aqcd 

While the benchmark model may appear to be a very specific choice for the hypercolor 
sector, its collider phenomenology is actually representative of a wide range of choices. In 
fact, the main effect of the choice of the hypercolor group and the representation to which 
the ip belong (encoded in the effective parameters in eq. (2.2)) can be represented by a small 
number of independent quantities, namely the masses of the coloron and the hyperpions, the 



^For a first attempt to simulate SU{N) gauge theories for arbitrary fcrmion representations and generic 
numbers of colors N, and in particular SU{2) with adjoint fermions, see ref. [9] and references therein. 
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production cross section of the coloron and its branching fraction into dijets vs. four jets (via 
hyperpions). It is therefore useful to describe the coUider phenomenology in terms of these 
quantities and we will now elaborate on how the phenomenology is affected as they are varied 
away from the benchmark model. For values of — ^ or — that are too large, the hyperpion 
decay mode may be suppressed or kincmatically inaccessible and dijets will be the dominant 
decay mode, on the other hand if the value of ^ is very small, then the hyperpions from the 
coloron decay will be very boosted and this can lead to the two jets from each hyperpion to be 
reconstructed as a single jet, thereby also turning the phenomenology effectively into the case 
where dijets are the dominant decay mode. It was shown in ref. [1] that the main constraint 
on the models at hand comes from resonance searches in the dijet channel, and colorons 
which decay dominantly into dijets are excluded at sub-TeV scales [10, 11, 12].^ On the other 
hand, as long as one stays away from these extreme regions the model is consistent with 
all present constraints and the collider phenomenology is virtually unchanged. Furthermore, 
above rrip ~ 1 TeV even a coloron decaying predominantly into dijets is not excluded by 
Tevatron dijet resonance searches, however we will not discuss this possibility in detail since 
there arc already dedicated searches for dijet resonances proposed for the LHC [13, 14]. We 
will use the parameters of the benchmark model to generate Monte Carlo signal events in 
section 3 simply because this is a representative choice, however our analysis strategy will be 
blind, and free of any assumptions regarding the values of the effective parameters. 

As mentioned, the value of the parameter % cannot be obtained from the QCD analogy 
(this term vanishes for QCD when G/^i, is replaced by F^^^), however NDA predicts its natural 
size to be 0(1). In fact, we can cross check this by considering similar couplings that appear 
in other models of new physics containing a massive spin-1 color octet state, such as a two-site 
model where the SU{3)l x SU{3)ji is broken by the vev of a scalar bifundamental and SM 
color is the unbroken SU (3) (The heavy vector adjoint in such models is known as an axigluon 
[15]) or the case of five dimensional QCD where the extra dimension has been compactified 
on a circle. In both cases the value of x is fixed at x = 1- In the case of 5-D QCD which is 
perturbative in the UV, this can be traced back to 5-D Lorentz invariance, and in fact x = 1 is 
the unique choice for which the scattering amplitude for G^^'^G^^^ G^^'^G^^^ remains unitary 
at high energies. As we will see in section 3, the main effect of x 'wiU be on the coloron pair 
production cross section, and we will adopt the choice % = 1 for generating signal events. 

The parameter ^ is another undetermined number of order one, and while we will demon- 
strate how the cross section for hyperpion pair production depends on ^, we will stick to the 
choice of ^ = in generating signal events, as this is the most conservative choice (non- 
zero values of ^ enhance the cross section). In terms of coloron decay branching fractions, a 
nonzero value of ^ gives rise to a two-gluon decay mode, but the branching fraction is quite 
small, about 7.5 x 10~^ for ^ = 1 so the collider phenomenology is unaffected unless ^ is much 
larger than its natural size. 

*A range of other potential constraints on this type of new physics were discussed in ref. [1] and shown to 
be irrelevant even at sub-TcV scales. Wc will not repeat this analysis here as we will be interested in TeV-scale 
masses where the constraints are even weaker. 
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Figure 1: Hypcrpion pair production cross section at the LHC, assuming the coloron mass is given 
through eq. (2.5). The green curve represents the resonant coloron production cross section which 
subsequently decays to a pair of hyperpions, while the red curve represents the pair production cross 
section of hyperpions from two initial gluons. The blue and violet curves demonstrate the enhancement 
of the gluon initial state cross section as ^ is allowed to be as large as its natural size. 

3. Search Channels and Strategies 

The most obvious place to search for the new states is the resonant production of the col- 
oron. This channel has strong discovery potential at the Tevatron for masses accessible there 
[1], however there are qualitative differences for the LHC. The most straightforward way to 
resonantly produce a coloron is through a q-q initial state, however the PDF's at the LHC 
suppress the resonant cross section in comparison to the Tevatron. While there is also the 
potential of resonant coloron production from a two gluon initial state, the cross section for 
hyperpion pair production is actually dominated by the minimal QCD coupling of hyperpions 
to gluons. This is displayed in figure 1, where the green curve represents the hyperpion pair 
production cross section coming from the coloron resonance, while the red curve represents 
the hyperpion pair production cross section from initial state gluons, with = 0. The blue 
and violet curves demonstrate the effect of turning on a nonzero value for ^, which is to in- 
crease the cross section. The most important consequence is that even if we can isolate signal 
events by the fact that they contain two resonances with equal masses, in general we will not 
be able to reconstruct a primary resonance. Thus at the LHC the coloron is more difficult 



Figure 2: A visual representation of the search strategy we propose. While the extraction of the 
hx peak may become unrealistic if the systematic errors of our background estimates are C(10), the 
success of the rest of the analysis should be robust. 

to discover compared to the hyperpions. After describing in section 3.1 the tools used in our 
analysis, we will demonstrate in section 3.2 the discovery potential for hyperpions in the four 
jet channel, then in section 3.3 we will augment our analysis by a search for pair produced 
colorons in an eight jet channel. We will demonstrate the range of masses where the LHC 
has a strong discovery potential by working with two cases, rup = 750 GeV (below which the 
coloron is discoverable at the Tevatron [1]) and nip = 1.5 TeV (above which in the benchmark 
model the coloron pair production cross section becomes too small, even though there is still 
significant room for the discovery of the hyperpions). 

Monte Carlo based methods for background estimates are expected to suffer from sys- 
tematic errors at high parton multiplicity, and we will take care in making statements about 
the eight-jet channel conservatively, using only cross section information for background (with 
expected systematic errors of order one) using the most sophisticated Monte Carlo tools avail- 
able for high multiplicity. We will argue that the large signal cross section will allow us to 
draw conclusions that should remain valid even in the presence of sizeable corrections to our 
background estimates. We will demonstrate in appendix A that our conclusions are indeed 
valid when applied to a plausibly realistic sample of background events, generated with our 
limitations of available computing power. 

Section 3.4 will be devoted to showing how this type of signal can be self calibrating 
and proposing a search strategy where one can extract the optimized cuts for the eight jet 
analysis using the mass peaks in the four jet analysis and the data at high hx- This procedure 
is summarized visually in figure 2, while we will argue the success of the four and eight jet 
analyses are robust, the extraction of the Ht peak may break down for systematic errors 
of 0(10) on the magnitude of our background estimates. The second step however should 
be regarded as optional and should not impact the discoverability of the hyperpions or the 
coloron. 

3.1 Simulation of Signals and Backgrounds 

Before discussing in detail the individual search channels and establishing the methods to 
extract signal from background we want to describe in detail the tools we used to simulate 
the signal processes and QCD multijet production. 
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3.1.1 Signal simulation 

For the calculation of signal cross sections and the generation of signal events at parton 
level we have incorporated the phenomenological Lagrangian presented in eq. (2.2) in the 
Sherpa Monte Carlo program [16] by implementing all interaction vertices of the model into 
Sherpa's matrix element generator AMEGIC++ [17], thereby allowing us to generate arbitrary 
tree-level processes within this model. We limit the use of full matrix-element calculations to 
the point of production of the hyperpions thereby capturing non-resonant QCD contributions 
as well as finite- width effects for the intermediate heavy states. The latter is of particular 
importance when studying coloron production as for our benchmark model the vector states 
are rather broad. In fact the ratio olTp/rrip ^ Q.2 prohibits the use of a simple narrow- width- 
approximation (NWA) for the coloron decays. Furthermore, using exact matrix element 
calculations, spin-correlations between the coloron decay products are properly taken into 
account. For the decays of the hyperpions themselves we use a NWA (r^f/m^j- 2 x 10^'^) 
and distribute the decay gluons isotropically in phase space consistent with the decay of a 
scalar intermediate state. The parton-level computations are done using the CTEQ6L parton- 
distribution functions [18]. The factorization and renormalization scales entering are set equal 
to the mass of the pair-produced states, namely hf = l^'R = m\ for hyperpion pair production 
and = jJ-R = fn^p for the coloron-pair channel. To account for additional jet radiation and 
intra-jet evolution the parton-level events are passed to PythiA (version 6.4) [19] for parton 
showering and subsequent hadronization. Hadron-level events are processed through PGS 
(version 4.0) [20] to include effects of finite detector resolution as well as jet reconstruction. 
We use a cone jet algorithm with i? = 0.5 throughout our analysis. 

3.1.2 Background simulation 

Simulating the Standard Model backgrounds for hyperpion and coloron pair-production im- 
plies generating inclusive four- and eight jet production, respectively. For the signal processes 
we expect rather hard jets from the decays of the heavy resonances, therefore we will always 
use cuts to demand high pT jets. Realistically estimating the rates for producing several 
well-separated and hard jets requires the use of a complete matrix-element calculation rather 
than the parton-shower approach. Lacking full one-loop QCD results we have to rely on 
tree-level calculations of the corresponding amplitudes. Intrinsically these computations suf- 
fer from significant scale uncertainties that affect both total cross sections and kinematical 
distributions. However, in our four jet analysis we will demonstrate very clearly that the 
hyperpion can be reconstructed such that the resonance is unmistakable even allowing for 
order one errors on the Monte Carlo background estimates. Coloron pair production, which 
is a 2 2 process, has a large enough cross section to overcome the 2 — >^ 8 background even 
without taking advantage of detailed kinematic features of the signal. We will further show 
that signal events can be reconstructed to pass very nontrivial kinematic cuts which further 
reduces background by a very significant amount. 

To compute the cross section of the four- and eight jet channels the new matrix-element 
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aAj{pT,j > 100 GeV) [pb] 


aij{pT,j > 200 GeV) [pb] 


^{ph) 


1463(13) 
2404(24) 
995(9) 


37.1(3) 
58.5(6) 
24.7(2) 



Table 1: LHC four jet production cross sections for jet-transverse-momentum thresholds of 100 

GeV and 200 GcV. The central factorization and renormalization scales arc given by the averaged 
transverse-momentum squared of the final-state jets. The phase-space cuts given in eq. (3.1) have 
been imposed. The numbers in parentheses indicate the absolute statistical error on the last digit. 

generator COMIX [21], interfaced with the Sherpa event-generation framework, is used. 

Comix constructs the QCD amplitudes using color- dressed Berends-Giele recursion [22] 
rather than traditional Feynman diagrams. To integrate the multi-parton phase space COMIX 
builds suitable phase-space mappings that are managed using an adaptive multi-channel 
method [21]. The combination of these techniques makes COMIX very efficient and fast and 
therefore well suited for our purposes. For all parton-level background calculations we use a 
dynamical factorization and renormalization scale, namely the average transverse-momentum 
squared of all final-state partons. In Tab. 1 we present four jet cross sections for two jet- 
Pt thresholds, namely pxj > 100 GeV and pxj > 200 GeV. Furthermore, all partons are 
required to satisfy 

ARjj>0.5 and |r/j|<2.0. (3.1) 

To estimate the scale uncertainty of the cross sections we consider variations of the factoriza- 
tion and renormalization scales by common factors of 1/4 and 4, corresponding numbers are 
also given in Tab. 1. The scale dependence of the cross sections can be estimated to about 
±50% and is dominated by the high power of the strong coupling entering the calculations. 
All results have been con&med with Alpgen [23]^. 

When considering eight jet production there is an enormous number of partonic subpro- 
cesses that contribute, however, most of them contain multiple pairs of quarks. Compared to 
the purely gluonic process these configurations are color-suppressed. To simplify our back- 
ground calculations we only include processes with at most two quarks in the initial and final 
state — all the remaining partons being gluons. The complete list of channels we take into 
account is 

gg^^g, gg^ 6g2q , gq 7glq , qq ^ 8g , qq Qg2q , (3.2) 

where q denotes all combinations of light-flavor quarks and anti-quarks allowed by QCD. As 
explicitly shown in ref. [21] this captures the bulk of the total eight jet cross section and we 
estimate the contribution of the neglected channels to be less than 20% of the total cross 
section. In Tab. 2 we list corresponding results for eight jet production at the LHC both 

^However, in Alpgen not all partonic subprocesses contributing to four jet production are available. There- 
fore a tuned comparison has been performed and processes with more than two quarks in the final state were 
excluded in COMiX as well. The contribution of these color-suppressed channels is less than 10% of the total 
four jet cross sections. 
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asjipTj > 100 GeV) [fb] 


asjipTj > 200 GeV) [fb] 


liPh) 


112(1) 
267(3) 
53.4(7) 


0.369(5) 
0.848(9) 
0.174(2) 



Table 2: LHC eight jet production cross sections for pxj thresholds of 100 GeV and 200 GeV. Note 
that only the constrained set of processes Usted in eq. (3.2) is considered. The appUed phase-space 
cuts are given in eq. (3.1). 

for our default scale choice and variations of that by factors of 1/4 and 4. The phase-space 
cuts used for this channel are again given by eq. (3.1), supplemented with pr thresholds 
PT,j > 100 GeV and pxj > 200 GeV. For these particular setups the scale variations change 
the results by approximately a factor of two up and down. This can largely be ascribed to 
the renormalization scale dependence of the 0(a|) process. 

For both signal and background, the parton-level events used in the four jet analysis have 
been interfaced to Pythia and PGS to account for parton showering and hadronization, as 
well as detector effects. In the eight jet case we only use total integrated cross sections to 
quantify background (except for the hx distributions presented in section 3.4) as we wish 
to avoid relying on detailed kinematic features of Monte Garlo based eight jet background 
events. Furthermore, the unweighting efficiency for eight parton phase space is rather low 
which makes the generation of large enough background event samples computationally very 
expensive. Once we demonstrate that the eight jet cross section for signal is larger than that 
of background at parton level after using simple pr cuts, we will proceed with a signal only 
study for mass reconstruction (after passing signal events through Pythia and PGS), which is 
justified since the small expected efficiency for background events to pass our reconstruction 
cuts should provide a nearly pure signal sample. To corroborate our statements about the 
dctcctability of coloron pair production in the eight jet channel we will confront the signal 
with a sample of showered and hadronized 6-parton background events. The corresponding 
analysis can be found in appendix A. 

3.2 Four Jet Analysis 

In order to improve signal over background, we take advantage of the main kinematic differ- 
ence between the two, namely that the jets in the signal events come from the decay of two 
equal-mass resonances and are therefore likely to have comparable and large pt's while the 
jets in the background are produced from QCD where the bulk of the background exhibits a 
strong pt hierarchy. Therefore for event selection, besides the generic cuts given in eq. (3.1), 
we impose a hard pT cut on the hardest four jets in each event. 

We then pair the four hardest jets in all possible ways and discard the event if there is 
no pairing in which the invariant masses of the two pairs are within 50 GeV of each other. 
If there is more than one such pairing, we take the one that yields the closest masses for the 
pairs. We then take the mean value of the pair masses as the candidate hyperpion mass. 
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and we also plot the correlation between the pair mass and the total invariant mass of the 
four jets. This separates signal from the bulk of the background because of the following 
reason: The background events that yield comparable pair masses have a special kinematic 
configuration, namely the jets that make up the pair that does not contain the hardest jet are 
quite far apart in AR in order to be able to match the invariant mass of the pair containing the 
hardest jet. This means that for a given value of the four jet invariant mass, the background 
events typically have a larger value for the pair invariant mass compared to signal events. 
This can be seen clearly in figures 3 and 4. 

3.2.1 Lighter mass case 

For nip = 750 GeV (and = 225 GeV) we choose a pT cut of 150 GcV on the four hardest 
jets. After the pairing and invariant mass cuts, the signal cross section is 2.2 pb and the 
background cross section is 24 pb. In figure 3 one can see that even though the pair mass 
has a peak at the hyperpion mass, the four jet mass values do not accumulate at the coloron 
mass, that is because in most of the signal events hyperpion pairs are created through their 
gauge couplings to gluons. 

We define the statistical significance in a histograms as 



where the bin size is chosen such that each bin is reasonably well populated while any 
significant shapes in the distributions are retained. We check that Xsig is invariant under 
changes in the bin size. For 1 fb""^ of data the statistical significance of the histogram 
in figure 3 is Xsig ~ 38. 

Figure 3 should also alleviate most concerns about systematic errors in using Monte Carlo 
backgrounds. As the pT cut that goes into this analysis is changed, the effect on the signal 
and background are quite different. The background distribution, in particular the position 
of the artificial peak introduced by the pT cuts, shifts as the cut value is changed. The signal 
on the other hand is very clearly peaked at the true mass of the hyperpion and only the height 
of the peak is affected as the px cut is varied, since the number of events passing the cuts 
changes. Therefore even allowing for errors of order one on the background, there will be a 
wide range for the value of the pT cut where the narrow resonant shape of the signal will be 
clearly distinguishable from the background. 

3.2.2 Heavier mass case 

For rup = 1.5 TeV (and mfr = 450 GeV) we impose a pr cut of 250 GeV on the four hardest 
jets. After the px cuts and the pair invariant mass cut the cross section of signal is 0.15 pb 
while the background cross section is 1.0 pb. In figure 4 we show that the signal still has 
very strong significance for the heavy mass case, for 10 fb^^ of integrated luminosity we get 
Xsig ~ 45. We reiterate at this point that there is strong discovery potential for hyperpions 




(3.3) 
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Figure 3: Upper plot: Pair invariant mass vs. four jet invariant mass for 1 fb~^ of data, red is signal 
while blue is background. The coloron mass is 750 GeV and the hyperpion mass is 225 GeV. Lower 
plot: The projection onto the pair invariant mass. 

heavier than 450 GeV as well, however in the benchmark model beyond this point the coloron 
pair production cross section is small and the analysis described in the next section becomes 
unrealistic. 

3.3 Eight Jet Analysis 

Even though we were able to demonstrate strong discovery potential for the hyperpions, we 
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Figure 4: Upper plot: Pair invariant mass vs. four jet invariant mass for 10 fb~^ of data, red is signal 
while blue is background. The coloron mass is 1.5 TeV and the hyperpion mass is 450 GeV. Lower 
plot: The projection onto the pair invariant mass. 

have not yet conclusively presented evidence for the existence of the coloron itself. For this 
we search for coloron pair production. The pair production cross section is plotted in figure 5. 
Note that the parameter x contributes to both the three-point as well as four-point couplings 
of the coloron to gluons, thus the coloron pair production amplitude has a piece proportional 
to s-iid the cross section can grow as x^- As discussed in section 2, we adopt the choice 
X = 1 as it appears to be the most plausible one in regard of similar theories containing 
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Figure 5: Coloron pair production cross section as a function of nip. The various curves represent 
the impact of different choices of the parameter x- 

color octet spin-1 particles. We can see that the cross section is quite large which allows 
us to compete with QCD backgrounds in terms of total cross section even before we impose 
kinematic cuts aimed at reconstructing the kinematic information in the signal events. 

We will present our analysis of the eight jet channel relying only on parton level cross 
section information for background, using a full 8-parton final state integration obtained with 
Comix. Generating a statistically large enough background sample of unweighted 8-parton 
events is unfortunately beyond our computational means, we will however confirm the validity 
of our conclusions in appendix A using inclusive 6-parton events that were then showered using 
Pythia, keeping in mind of course that kinematic shapes in the background are expected to 
have significant systematic uncertainties. 

3.3.1 Lighter mass case 

In order to decrease the background cross sections as much as possible we use an incremental 
Pt cut scheme. We determine the cut values from a Monte Carlo study of signal events, 
we choose a px cut for the eight hardest jets such that the efficiency of the cut on the 
eighth hardest jet is ~ 55% and each additional cut has a relative efficiency of ~ 80%, this 
sets the efficiency of the combined cut scheme on the signal to be ~ 13%. In particular for 
mp = 750 GeV we demand pTjs > 40 GeV, prji > 60 GeV, prje > 90 CeV, pT,j5 > 125 GeV, 
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PT,j4: > 160 GeV, pT,j3 > 200 GeV, pT,j2 > 250 GeV, prji > 320 GeV. With these cuts the 
background cross section is 1.2 pb at parton level, and is expected to decrease further once the 
cuts are imposed at jet level, since parton showers will typically split the jets and therefore 
reduce the pt of individual jets (this expectation will be confirmed in appendix A). The 
signal cross section after imposing the pr cuts past parton showering and jet reconstruction 
is 2.6 pb. 

Despite the systematic uncertainties in computing the eight jet backgrounds, these num- 
bers indicate that even without taking advantage of any kinematic information in the events 
such as invariant masses of pairs of jets, signal and background cross sections are expected to 
be comparable in size even when one allows for relatively large if-factors. Below we present 
a signal only study to demonstrate that with the proper analysis cuts we can reconstruct the 
cascade decays in the signal events, while the efficiency for background events to pass these 
reconstruction cuts is expected to be very low. This will be confirmed for showered 6-parton 
events in appendix A. 

First we do a Monte Carlo study of signal events in order to compare the parton level 
truth information to the result of parton showering, hadronization and detector effects. Wc 
study the fate of the eight jets coming from the decay of the coloron pair and we find that 
fairly often parton showering will split the jets and make the original kinematics very hard 
to reconstruct. Furthermore, with the high number of jets present in the final state it is 
not uncommon for at least two of the primary jets to partially overlap and be reconstructed 
as a single jet. We do find however, that in the events where the eight primary gluons are 
faithfully found as eight separate jets, they are virtually always among the ten hardest jets 
in the event, and the hardest four jets almost never come from additional radiated partons. 
Based on this information, we design our analysis procedure as follows: We take the ten 
hardest jets in an event (if there are only eight or nine jets in the event, we take those, while 
the rest of the procedure is unchanged), we then take all possible subsets of eight jets out of 
those always keeping the hardest four, and wc pair them up into four pairs in all possible ways. 
We then demand that all four pairs have an invariant mass in a window that corresponds to 
the signal peak in the four jet analysis (the histogram in figure 3). In the lighter coloron case 
this is 175 GeV < m2j < 245 GeV. We drop events for which no such pairing exists. We 
expect this cut to severely reduce background as background events normally display a large 
hierarchy of pt's and it is very unlikely that four pairs with similar masses can be formed. 
In fact, considering the small pairing efficiency of background in the four jet analysis where 
wc demand only two pairs with similar invariant masses, we are confident that requiring four 
pairs with comparable masses will yield a virtually pure signal sample. For the background 
analysis on 6-parton events presented in appendix A this is shown to be valid. 

The signal cross section after the invariant mass cuts is found to be 0.47 pb. While most 
events that pass the cuts have a unique pairing that satisfies the cuts, there are events with 
more than one possible pairing, the distribution is given in table 3. Instead of throwing away 
the events with more than one possible pairing we choose to assign a weight to all possible 
pairings in an event that pass the invariant mass cuts (i.e. where all four pairs fall in the 
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cr with 


1 pairing 


2 pairings 


3 pairings 


more pairings 


(in pb) 


0.28 


0.10 


0.04 


0.05 



Table 3: Cross section (after px cuts) of signal events with one or multiple candidates for four 
dijet-pairs all compatible with 175 GeV < m2j < 245 GeV. 

correct mass window) and take a weighted average for the coloron masses as follows: For 
any pairing V that passes the invariant mass cuts (as an example V could correspond to jets 
(1,5), (2, 9), (3, 4), (6, 8) being paired up), let {mi^-p^i} be the values of the four reconstructed 
hyperpion masses. Now as we pair the hyperpions themselves, we get three possible pairs of 
coloron masses, let these be ("T'p<,(p,j) > "^p>,(P,j))5 where j = 1,2,3 (continuing with the above 
example, the pair (1,5) could be paired with (2, 9), (3,4) or (6,8)). We define the (unnormalized) 
weight for any such pairing as 

^{V,3) = 1 > (3-4) 



where 



1 



1 



i=l 



On 



P,{r,j) 



2 



\^p,ir,j)) 



(3.5) 



i=l 



This weighting scheme favors pairings for which the hyperpion masses come out close to each 
other, and similarly it favors the pairing of the hyperpions which yields coloron masses that 
are close to each other. 

For each event we then calculate one set of coloron masses as the weighted average 



y^p>l — 



(3.6) 



The results are plotted in figure 6 for 5 fb~^ of data. The bulk of the signal is indeed 
reconstructed near the true mass for the coloron. Note that the main input for making this 
plot is the extraction of the hyperpion mass from the 4-jet analysis, a feature of the signal. 
Thus there is no reason to expect the background to have an accumulation point, the showered 
6-parton background events in appendix A confirm this expectation. 

3.3.2 Heavier mass case 

We use the same techniques for the case of the heavier coloron. Even though we are pair 
producing 1.5 TeV colorons, we find that the signal cross section is large enough for this 
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Figure 6: Coloron pair masses calculated from signal events with (5 fb ^) of luminosity accumulate 
near the true mass. The two circles represent IFp and 2Tp distances from the true mass. 



a with 


1 pairing 


2 pairings 


3 pairings 


more pairings 


(in fb) 


1.1 


0.32 


0.14 


0.17 



Table 4: Cross section (after px cuts) of signal events with one or multiple candidates for four 
dijet-pairs all compatible with 350 GeV < m2j < 475 GeV. 

analysis to be possible. The pT cut scheme we are using for the heavier mass coloron is given 
by PT,js > 75 GeV, prji > 115 GeV, prje > 180 GeV, pT,j5 > 240 GeV, ptja > 300 GeV, 
PT,j3 > 380 GeV, PT,j2 > 470 GeV, prji > 600 GeV. With these cuts, the background cross 
section at parton level is 5.7 fb, which once again is expected to drop further once the cuts 
are imposed after parton showering and jet reconstruction. The signal cross section after 
showering and jet reconstruction is 9.5 fb. 

We then use the mass window 350 GeV < m2j < 475 GeV (in accordance with figure 4) 
to veto events in which the jets do not pair up correctly, as described for the lighter mass case. 
Once again, the majority of the events past the invariant mass cuts have a unique pairing, as 
demonstrated in table 4. Using the same technique as in the lighter mass case, we find that 
with 200 fb"^ of data one can reconstruct signal events quite successfully. Our results are 
plotted in figure 7. 
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Figure 7: Coloron pair masses calculated from signal events (200 fb ^) accumulate near the true 
mass. The two circles represent IFp and 2Tp distances from the true mass. 

3.4 Guidelines for a Self Calibrating Search 

While in sections 3.2 and 3.3 we have highhghted the discovery potential of the LHC for the 
hyperpions and colorons for a wide range of masses, the px cuts in each case were chosen 
to enhance signal over background. In the following, we will attempt to prescribe a search 
strategy that can be used in an experimental search without prior knowledge of the scale of 
new physics. We will also address which stages of the analysis could be susceptible to large 
systematic errors in our background estimates relying on tree level amplitudes. 

The only input entering the four jet analysis we presented in section 3.2 was the pT cuts 
on the four hardest jets. (We find that demanding the two hyperpion candidate masses to 
be within 50 GeV of each other works for the entire mass range of interest.) However, the 
prominence of the hyperpion mass peak above background in figures 3 and 4 demonstrates 
that this analysis can be done with a sliding pT cut. As the cut value is increased, the 
background distribution will shift to higher energies and the signal peak separates cleanly 
from the background for reasons described in section 3.2. Thus the four jet analysis can be 
modularized, and the position of the hyperpion mass peak can be used as an input to the 
eight jet analysis, where we demand four pairs of jets with an invariant mass in the window 
determined by the four jet analysis. Note also that this procedure is robust to errors of 0(10) 
in the background estimates, since Xsig scales as therefore even with a i^T-factor of 10 
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Figure 8: The hx distribution of signal (red) and background (blue) with a cut on the px of the 
eighth hardest jet. Signal cross sections represent coloron pair production with mp — 750 GeV. From 
left to right, the value of the pt cut is increased from 50 GeV to 75 GeV to 100 GeV. 




Figure 9: The hr distribution of signal (red) and background (blue) with a cut on the px of the 
eighth hardest jet. Signal cross sections represent coloron pair production with nip = 1.5 TeV. From 
left to right, the value of the px cut is increased from 100 GeV to 150 GeV to 200 GeV. 

the statistical significance in both our benchmark studies remains above 10a. 

Note that this way of performing the four jet analysis is very general and does not assume 
the existence of the coloron. Since the signal cross section is dominated by the contribution 
from minimal QCD couplings of a color octet scalar, we have demonstrated strong discovery 
potential for any state which is pair produced with a comparable cross section and which 
dominantly decays to dijets, with the implicit assumption that the coupling giving rise to 
this decay mode is small, otherwise this state may first be observed as a dijet resonance. A 
similar analysis for states appearing in non-minimal technicolor models has been performed 
in ref. [24]. Spin-1 states decaying to dijets and scalars decaying to heavy flavors have been 
looked at in ref. [25] . Even though ref. [26] investigated pair production of pseudoscalars that 
predominantly decay to four jets, they concentrated their efforts on a subdominant decay 
channel including a photon. The phenomenology of ref. [27] also includes a resonance channel 
with a four jet final state, however they focus in their analysis on final states with top quarks. 

To demonstrate how the optimized pT cuts entering the eight jet analysis can be obtained, 
we focus on the hx distribution of eight jet events with a sliding px cut on the eighth hardest 
jet. (/ly is defined as the scalar sum of all jet transverse momenta.) The results are displayed 
for the lighter and heavier coloron cases in figures 8 and 9. At a particular value of the pT 
cut we find that signal becomes visible as a bump above background (left plot) with a peak 
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Figure 10: The peak position of the Ht distribution in the signal as a function of the coloron mass, for 
different values of the px cut on the eight hardest jet. The behavior is linear to a good approximation. 
The error bars indicate the uncertainty in specifying a precise position for where the distribution 
peaks, and become larger as the number of events in the distribution decrease as a consequence of the 
increasing px cut. 



that is shifted away from the artificial peak in the background introduced by the cuts. As 
the pt cut is increased further, signal falls less sharply than background (middle plot), and 
while at large values of the cut signal overwhelms background in cross section, the difference 
between signal and background becomes one of normalization rather than one of shape (right 
plot) since the pT cut begins to shift the signal peak as well. While these plots are made 
at parton level, we expect parton showering to affect signal and background shapes similarly 
since there should be no significant difference in the way they radiate, and our conclusions 
should still apply. The point we wish to make here is that even if the background estimates 
needed to be scaled up by a i^-factor of 0(1), the signal would still be visible as a bump at 
a higher value of the pT cut (e.g. in the middle plot rather than the left plot) due to the 
different rate at which the signal and background cross sections fall with the increasing cut 
value. This procedure can break down if the i^T-factor is so large (0(10)) that the signal cross 
section would only become comparable to background by the time that the pT cut is large 
enough to shift the signal peak position to match the background shape (right plot), in which 
case there may be no clear indication of a signal bump. 
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Figure 11: A scheme for picking the incremental pT cuts on the eight hardest jets as a function of 
the coloron mass. The values are chosen such that the efficiency of the cut on the eighth hardest jet on 
signal events is 55%, and the efficiency of each consecutive cut is 80%, such that the total efficiency of 
this cut scheme on signal is 12%. We find that this choice gives a higher S/B enhancement compared 
to a single pt cut on the eight hardest jets. 

Assuming that the systematic errors in our background estimates are not so large as to 
invalidate the procedure described above, the position of the hx peak can be used to deduce 
the coloron mass, as they are linearly correlated. This is displayed in figure 10 where we plot 
the peak position of the hx distribution in the signal as a function of the coloron mass using 
events that have been passed through Pythia and PGS. While this plot has a dependence 
on the jet reconstruction algorithm that is used as well as details of detector performance, it 
demonstrates that a linear relation exists between the signal peak of the Ht distribution and 
the mass of the coloron. We expect that any experimental search would reproduce figure 10 
using their own jet reconstruction algorithm and detector simulation. 

This procedure makes it possible to link the excess in the hx distribution to a candidate 
mass for the coloron and use that to optimize the incremental pT cuts that enter the eight 
jet analysis. Note that our analysis does not assume ^ = 0.3 and is therefore applicable 
to more general scenarios than the benchmark model. In figure 11 we display a scheme for 
choosing the incremental pT cuts that enter the eight jet analysis, the values of the px cuts are 
chosen for each value ol mp to have the same efficiency on the signal as described in section 



3.3. Combining this with the result of the four jet analysis to set the pairing mass window, 
the eight jet analysis can now be performed. 

Note that while we have proposed using the hr distribution to determine the px cuts to 
be used in the coloron mass reconstruction analysis, this is not an essential step. Even with 
large errors on the background, we have shown that the eight jet signal and background cross 
sections are comparable before any reconstruction cuts are imposed, which are expected to 
further enhance signal over background significantly, this is confirmed in appendix A with a 
background sample of showered 6-parton events. With this in mind, it should be possible to 
perform the coloron reconstruction analysis with a sliding px cut on the eight hardest jets 
instead of our optimized incremental pT cut scheme. While this may decrease S/B by an 
0{1) factor, the overall success of the analysis presented in section 3.3 should be robust. 

4. Conclusions 

Even though the production of colored states is the main strength of hadron colliders, at the 
present there exist no dedicated searches for new physics in multijet channels because of the 
large background cross sections and the presence of systematic uncertainties in background 
estimates based on tree level calculations. Continuing in the spirit of ref. [1] , we have reiterated 
in this paper that there is a family of theoretically very simple models, relying on a minimal 
set of assumptions on how the SM is extended, with a subtle phenomenology where the 
states of new physics appear almost exclusively in multijet channels. While these models can 
incorporate diverse choices for the new physics sector, the existence of colored octet vectors 
and scalars as bound states of the underlying dynamics is generic. 

We used a phenomenological Lagrangian to describe the production and decay of these 
states. There is a choice for the microscopic physics for which the strong dynamics are an 
exact copy of QCD which allows us to extract the effective parameters of the phenomenolog- 
ical Lagrangian. We utilized this benchmark model to show quantitatively that the collider 
phenomenology is interesting and argued that this conclusion remains generic for variations 
around this choice. We then demonstrated for two choices of the coloron mass that hyperpi- 
ons can be discovered in the four jet channel with relatively low integrated luminosity. We 
have furthermore shown that colorons can be discovered in the eight jet channel with ad- 
ditional statistics. While we lack the computational power to simulate a statistically large 
enough sample of unweighted 8-parton background events, we did demonstrate for a sample 
of 6-parton events with subsequent showering that our expectations are confirmed. 

We have further argued that the signal can be self calibrating and we have proposed a 
realistic search strategy that assumes no prior knowledge of the scale of new physics, where the 
hyperpion search can be performed with a sliding px cut, and an excess in the hx distribution 
of eight jet events can be used to optimize the px cuts for coloron mass reconstruction. While 
this last step may be susceptible to large errors in our background estimates, it is not essential 
for the discovery of the coloron. 
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We remark that the colhder phenomenology studied here can also be used in a more 
agnostic context. One of the most obvious new physics scenarios to expect at a hadron 
collider is the pair production of a heavy colored state. Barring the case where this state 
is quasi-stable (which gives rise to R-hadron signatures) it can decay in a number of ways, 
and since the decaying particle is colored, the decay into dijcts is one of the most robust 
possibilities. Once again, with this in mind it is surprising that model-independent collider 
searches for new physics in multijets have not been performed in the past. We would like to 
emphasize that the four jet analysis we presented is independent from the existence of the 
coloron, thus we have in fact demonstrated the discovery potential for the pair production of 
any colored state which is produced with a cross section comparable with a color octet scalar 
and which decays to dijets, with the assumption that the coupling giving rise to this decay 
mode is small, otherwise the primary production mechanism would be resonant production, 
and this state would have first manifested itself as a dijet resonance. 

The LHC-period is a crucial time in particle physics, and while we may expect to discover 
a simple mechanism that solves the hierarchy problem, we should be prepared for discoveries 
in other channels which may (at first) appear unrelated to such a mechanism. The models 

presented here are simple, they can occur readily in gauge field theory, they are consistent 
with all experimental constraints and their phenomenology overlaps with the main strength 
of the LHC, the production of colored new particles. Whether or not they may be connected 
to a larger mechanism having to do with the electroweak symmetry breaking, they provide a 
template for a new physics scenario appearing in the multijet channel. In this work we have 
demonstrated that the LHC has strong discovery potential for such new physics, even in the 
case that it may arise from a different underlying mechanism than the one we have used to 
model our phenomenological Lagrangian on. Therefore we propose that multijet resonance 
searches, like dijet resonance searches, should become standard at hadron colliders. 
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cr with 


1 pairing 


2 pairings 


3 pairings 


more pairings 


(in fb) 


17 


6.8 


2.0 


2.6 



Table 5: Cross section (after px cuts) of background events with one or muhiple candidates for four 
dijet-pairs aU compatible with 175 GeV < m2j < 245 GeV. 



a with 


1 pairing 


2 pairings 


3 pairings 


more; pairings 


(in fb) 


0.14 


0.05 


0.01 


0.02 



Table 6: Cross section (after pt cuts) of background events with one or multiple candidates for four 
dijet-pairs all compatible with 350 GeV < m2j < 475 GeV. 

A. Check of Our Conclusions Using Background Events 

Even though we use COMIX for a fuU 8-parton phase space integration, generating a sample of 
unweighted 8-parton events is beyond our computational means. In order to test the vahdity 
of our conclusions we have therefore chosen to generate a sample of 6-parton events which 
were subsequently showered, hadronized and reconstructed using Pythia and PGS. Since 
these results should only be taken as a sanity check rather than a fully realistic background 
study, we choose to present them outside the main body of this paper. 

A.l Lighter Mass Case 

After the incremental pT cuts described in section 3.3.1, we find that the background cross 
section to be 1.0 pb, lower than the parton level result as expected. Upon demanding the 

presence of four pairs of jets in the mass window 175 GeV < m2j < 245 GeV we find that 
97% of the background is eliminated, the cross section of background events with one or more 
possible pairing is given in table 5. Using the reconstruction techniques described in section 
3.3.1 we superimpose the background events on top of the signal in figure 12. It is evident 
that background events do not have an accumulation point while the signal favors the true 
coloron mass. 

A. 2 Heavier Mass Case 

Using the incremental px cuts described of section 3.3.2, the background cross section is found 
to be 4.9 fb, lower than the parton level result in accordance with expectations. Requiring the 
presence of four pairs of jets in the mass window 350 GeV < m2j < 475 GeV eliminates 96% of 
the background, the remaining cross section being distributed as in table 6. The results after 
reconstructing the coloron mass are displayed in figure 13. Once again, signal accumulates 
near the true coloron mass while the background has no obvious point of accumulation. 
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